Entanglement dynamics under local Lindblad evolution 
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The phenomenon of entanglement sudden death (ESD) in finite dimensional composite open sys- 
tems is described here for both bi-partite as well as multipartite cases, where individual subsystems 
undergo Lindblad type heat bath evolution. ESD is found to be generic for non-zero temperature of 
the bath. At T = 0, one-sided action of the heat bath on pure entangled states of two qubits does 
not show ESD. 
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I. INTRODUCTION 

Entanglement is considered to be the most useful re- 
source in Quantum Information Theory [1] : it is essential 
for quantum teleportation, superdense coding, commu- 
nication complexity problem, one-way computation etc. 
Not only that creating entangled state is a non-trivial 
task, to store or transmit entangled states in an error- 
free manner is also difficult, if not impossible, due to the 
very fragile character of quantum system- every quan- 
tum system has a high possibility to interact with its 
environment, and thereby, the system will, in general, 
get entangled with its environment. This will give rise to 
the phenomenon of decoherence [2J. 

In general, the purity of any initial state of the quan- 
tum system goes down with time in the presence of deco- 
herence. This decoherence time depends on the system as 
well as on the character of the interaction of the system 
with its environment. So, due to the monogamy property 
of entanglement, the initial entanglement (if any) of a bi- 
partite or multipartite quantum system will, in general, 
decay (to zero) when each individual system undergoes a 
decoherence procedure. What can be said about the as- 
sociated rate of the above-mentioned decay in entangle- 
ment? How does one compare the rate of decoherence of 
the individual subsystems and the rate of decay in initial 
entanglement among the subsystems? In this connection, 
Yu and Eberly [3] described a phenomenon called "entan- 
glement sudden death (ESD)" in which the entanglement 
decay rate is shown to be exponentially larger than the 
rate of decoherence. This happens whenever the individ- 
ual qubits of a two-qubit system undergo evolution under 
local heat bath action at zero temperature. 

ESD was shown for a certain class of two-qubit states 
which were initially entangled at zero temperature (T = 
0). It was then generalized for T ^ for U X" states 
[4]. In this paper, we want to find out the largest set 
of two-qubit states each of which undergoes ESD due to 
Markovian local heat bath dynamics. If the initial bath 
state is the thermal state, at finite temperature (T ^= 
0) it is shown that all two-qubit states show ESD. The 
same is true also where initial states of the bath is a 



squeezed thermal state. On the other hand, under one 
sided Markovian heat bath evolution (i.e, only one of the 
two qubits is undergoing Markovian heat bath evolution) , 
one does not see any ESD at T = for any two-qubit pure 
state if the initial state of the bath is the vacuum state. 
However, we will find ESD if the initial state of the bath is 
the squeezed vacuum. This feature does not hold for one- 
sided quantum non demolition (QND) evolution |10| with 
squeezed thermal state as the initial state of the bath. In 
this case, no two-qubit pure entangled state show ESD. 
In deriving all this, we have used the factorization rule 
for concurrence [7]. We have also described the ESD 
phenomena for states of d<£> 2 systems as well as those of 
n-qubits by using the above factorization rule. 

The structure of the paper is as follows: in section [H] 
we discuss the ESD in two-qubit system with the bath 
acting only on one qubit. In section |III| we generalize 
the factorization law for entanglement decay [7] to d 2 
dimensional systems. In section |IV[ we use the result of 
section III and show ESD in n-qubit systems. In section 



|Vj we show ESD when the bath is in squeezed thermal 
state initially. In section |yi} we consider the QND-type 
evolution of the individual qubits. Finally in section |"VII| 
we will derive the sufficient condition for ESD in all di- 
mensions. Section IVIIII contains our conclusion and dis- 
cussions on some open problems. 



II. TWO-QUBIT UNDER LOCAL BATH 

We study a 2-qubit system initially prepared in an en- 
tangled state where one of them (say the first one) is 
interacting with a reservoir at temperature T. If the 
initial state of the bath in contact of the system under 
consideration be the thermal state, then the dynamics 
of a single-qubit density matrix p describing the system 
(under the Born-Markov-rotating wave approximation) 
is given by: 



dp_ _ {N+ 1) 7 
dt ~ 2 



[2<7_/9<7+ — (T + (T_p — p<7 + 0~-\ 
2cT + p<7- — CT-CT+p — pCT-(J + ] (1) 
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where N is the mean occupation number of the reservoir, 
7 is the spontaneous decay rate of the qubits, a+ = |1) (0| 
and <7_ = |0) (1|. Here we are ignoring the unitary part of 
the evolution which is irrelevant for our purpose [31 112) . 
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We can rewrite eq n ([TJ as: 

p = Hp] pij = 51 Lij,kiPki 



\x)ab{x 



(2) 



where L is the matrix representation for A (called Lind- 
blad operator [32] )■ The solution for eq"([2| is: 



(3) 



hi 



Here V = exp(Lt) [6] is a completely positive map as 
p(t) is a valid density matrix of the qubit at any time 
t. Our aim is to find entanglement in the evolved 2- 
qubit state after applying the map V on one qubit. 
For that purpose we will use the factorization law for 
entanglement decay [7] which states that the concur- 
rence of a two qubit pure state \x)- with one qubit be- 
ing subject to an arbitrary channel V— will satisfy the 
equation C(p AB ) = C (\ X )) C ((I ® F)(|0+}<^+|)) where 
Pab = (I ® ^0 (|x)(xl)- To show entanglement sudden 
death (ESD) in any state [3], it is sufficient to show that 
the state \cj> + ) evolves to a separable state under the ac- 
tion of given single qubit operation V. We can calculate 
the V matrix for our case: 



V 



a; cot (6>) 
x 

Ox 




yi 



!I2 

o 
o 

- x cot(6>) 



(4) 



where j/i = 2xcot(9)N, y 2 = 2x cot(6)(l+N), tan(6>) = 
( 2x{i+2N) ^ ^ x = exp [_i 7 (! + 2 JV)t] . When we apply 

this V on one side of |<^ + ) we will get the mixed state: 



M(t) = 



hxcot(0) 

2/2 

2/! 

x 



x 





1+x" 



xcot(6 l ) 



(5) 



For the matrix M(t) to represent a separable state, the 
partial transpose of M(t) should be a positive semi- 
definite matrix [BUS]. If sinh( 7 (l + 2N)t > 2^ ^+j] , 

the matrix M(t) is positive semi-definite under partial 
transposition and hence is separable. This shows that 
the operator V acting on one-qubit evolves the maxi- 
mally entangled state |^ + ) into a separable state. The 
factorization law for entanglement decay implies that all 
the two qubit pure entangled states show ESD. 

Now consider the zero temperature case (T — 0). For 
T = the mean occupation number N is zero and hence 
2/i is also zero. The matrix M matrix will never be pos- 
itive under partial transposition and hence will always 
represent an entangled state. The direct implication of 
this result is that there is no pure two qubit entangled 
state that will show ESD at T = 0. 
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FIG. 1: Duality between the pairs (|x)(xIj $) an< i ($x> P($)) 
via Choi-Jamiolkowski isomorphism. 



For T 7^ 0, every two-qubit pure state shows ESD 
and since mixed states are convex combinations of pure 
states, every two-qubit mixed state will show ESD. How- 
ever, at T = some mixed states do show ESD even 
though no pure state does. 

In section [TV] we will extend our result to the n-qubit 
case. But before that we will derive the generalized fac- 
torization law for entanglement decay, i.e, the factoriza- 
tion law for d®2 systems in the next section. 



III. FACTORIZATION LAW FOR 
ENTANGLEMENT DECAY FOR d <g> 2 SYSTEMS 



In this section, we consider the effect of any single- 
qubit trace-preserving map $ on side B of any bipartite 
density matrix pab, where dimHA — d. In ref. [T7] 
one can see some discussion on the factorization law for 
entanglement decay for d <E> d systems. 

To start with, we consider first its effect on any pure 
state \x)ab having Schmidt form \x)ab — y/p\00)AB + 
\/l — p\ll) ab- Here |0)b, |1)_b are the eigenvectors of ct 2 
for the subsystem B and 10)^, are taken from the 

standard orthonormal basis {|0)a, |1)a 3 1 ■ • , \d— 01 
T-La with < p < 1. 

As (1a <8> $) (|x)ab(xJ) is a two-qubit density matrix, 
one can find out its concurrence C((1a ® $) (Ix)ab(xI))- 
So, according to fig.[l], we have 



C((U®$)(|X>AB<Xl)) = 

C{M: (<P + \DE(\x)AD{x\®PEc($))\<i> + ) 



DE 



(6) 
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where p EC = <Je ® $) (|</>+>£c(0 + |) , M : Z := Z/(TrZ) 
and \4> + ) = 4|(|00) + |11)) is a two qubit maximally 
entangled state. 

We n ow perform the POVM {M = y/p\0) A (0\ + 
y/l-p\l) a (l\, I d xd - M} on side A of the d ® 2 den- 
sity matrix pac($) -So in the case when M is clicked, the 
normalized output state is given by 

M : (M®1 c )pac($) (M(g>l c ) ■ = oac say. (7) 

Note here that pac{§) is same as pec($) with just E 
replaced by A. Along the line of [TJ, it will follow imme- 
diately that 

N:(<h+\ 

DE (|x).4.d(xI ® PEC(§)) \ c I )+ )dE '■ — (8) 

Thus we see that (l A <g> $) (Ix)ab(xI) — &ac- So 



C{{I a ®$){\x)ab{x\)) = C{a AC ). Now the analysis of 
[7J will follow straight away and we will get: 

C((Ia ® $) (\x)ab(x\)) =C ((I a ® $)(|0+)ab(0 + |)) 

xC(\ X )ab(x\)- (9) 

Next we consider the case when |x)ab has the following 
Schmidt decomposition: 

\x)ab = (u A ® t/s)(Vp|oo>Ai3 + vT^piii)^) 

= ® ^b)|Xo)ab (10) 

where U A is a d x d unitary matrix while Ub is a 2 x 2 
unitary matrix. 



(I A ® $)(|x>abvX|) = -A/" : (M' ® 7 C ){(^ ® $)(I^>^(^I)}(X' ® /c) := (say) , (11) 
where .M' = U A MU\ = U A (y/p\0) A (0\ + V^Pl 1 ) aW) u a and ^t) ab = ( U A ® U B )\4>+) AB . One can then write 

g ac = M : (U A ® /s)[(M ® /o){(/A ® $)((Ja ® C^b )l</> + ) ab(4> + I (*A ® ^))}(M ® /c)](t/l ® Is) : (12) 



Now {I A ® ^b)|^ + ) as = (Va ® ^b)|0 + )^ b where the 
2x2 unitary matrix JTb is acting on the two dimen- 
sional space spanned by {|0) B , and the d X d uni- 
tary matrix Va acts on the linear span of {10)^,11)^} 
and is extended unitarily on the entire space spanned by 
{\0) A , \1) A , ...,\d- 1) A }. This means that both V A \0) A 
as well as Va|1)a lie inside the linear span of { 1 0) ^ , |1)aI- 
Note that the same is true for Ub also. 

So, here a AC = Af : {U A ® I B )[(MV A ® Ic){(Ia ® 
$)(\0 + )ab^ + \)}(v1M ® J C )](E^ ® J B ) := N : (U A ® 
/flhc(Pi ® Is) : 

Therefore, C(ctac) = C(t ac ). 

Now 



we have for any pure state |x) of a d £g> 2 system: 



C((Ia®$)(|x>ab(xD) 
=C ((Ia ® $)(|^+)ab(0 + |)) C(|xo)ab<Xo|) 
=C ((U ® $)(|^+)ab(0 + |)) C(|x)ab<xI) (14) 



Let pab be a mixed state on C d £§><C 2 . Let us consider 
an arbitrary ensemble representation for p A B- 



tacTac = T AC ((Jy ® o~y)T AC (a y ® <r y ) 
= (XVa ® Ic)/9ac($)(v1-M ® /c)K ® CT y) 

x (xy|®/ c )pAc($)*(vJx®/c)K®^), (13) 

which will, in turn, show that det [t A ct A C ~ 
A/ 4 x4] = det[r) A c{McryM ® <r y )r?i* AC (M<TyM ® <r y ) - 
A/ 4 x4] with ?7ac = (Va ® 7cWc(S)(Va ® 7c). 

Thus the concurrence of tac is same as that of the two- 
qubit state T] EC . So C(0U ® $) (\x)ab{x\)) = C(rj EC )- 
Now, for the state i]ec, we have already seen that the 
factorization rule holds good for concurrence. Therefore, 



JV 



Pab = y^ J Pj\(t>i)AB{4>j\, with (15) 
^f(pab) < 5>,£f(|^)ab(^|), (16) 



Ep(p) being the entanglement of formation. Now 



E F (\^) AB (^\) = /(C(I^-)ab(^D) , (17) 
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where f{x) is a monotonic function of x € [0, 1]. Here 
E F ((I A ®$)( PAB )) 

= E F (^PiOUSSXI^nfail) 



JV 



<X>^((^®$)(|&>Afl<&|)) 

J'=l 

JV 

= ^P,/(C((I A ®$)(|^) AB (^|))) 



JV 

= Y^Pjf (C (OU ® C (\4>j)ab(<Pj\) 

(18) 

Nevertheless, when C ((I a <8 $)(|0 + )ab(0 + |)) becomes 
zero / (C (QU <8 $)(|0+)ab(^ + |))C will au- 

tomatically become zero. So, in that case, Bp ((IL^ (g> 
$)(pab))=0. 



IV. ESD IN n-QUBIT SYSTEM 

In this section we will be dealing with the evolution 
of entanglement of an n-qubit system. We will prove 
that at finite temperature, all the n-qubit states show 
ESD under the map V- given in eq™|4])- acting on each 
qubit locally. To make the proof less cumbersome we will 
do it for the three-qubit case only. Let the state of the 
three-qubit system be \iP)abc evolving under the bath 
V <E) V <g> V. The final state is: 

P = {V®V®V){\iI))abc{4>\) 
= (I <g I ® V)(l <8 V <8 I)(V <8 1 <8 1) (I^Usc^l) 
= (I® I ® V)(I <8 V <8 1)pa-.bc 

= (I®I®V)pA:B:C (19) 

where pa-.bc = (V ® I ® I) (|^>abc(0|) and p^^c = 
(I (g) V <g> I)pA:BC- We have seen in section [TT] that V 
evolves to a separable state in some time t, and so 
the factorization law for entanglement decay for 2<gd sys- 
tem implies that the state Pa-.bc will be separable in the 
partition A : BC at time t. As this time t is the maximum 
time any state can take to lose entire entanglement by the 
action of channel V on one subsystem, we can see that 
at the same time t the reduced state pbc = Trf/j^Bc] 
becomes separable in the partition B : C. Hence we get 
the full separability in the partition A : B : C. So all 
the three-qubit pure states show ESD. Mixed states are 
the convex sum of pure states. If all the pure states show 
ESD, all the mixed states will also show ESD. This result 
can be generalized to any the n-qubit case and we will 
get full n separability at time t. 



V. SQUEEZED THERMAL BATH 

In this section we will consider the case of an n-qubit 
system S interacting with a squeezed thermal bath acting 
locally on each of the n individual qubits. The evolution 
of the reduced density matrix of the system S in the 
interaction picture has the form [T0l - fT2| : 



7 o(AT+ 1) U.p s (t)a + - ±a+a-p s (t) - \p s (t)a + o- 
+ l0 N (a + p s (t)a- - ^a-a + p s (t) - l -p s {t)<j_<j + 



-f Ma+p s (t)a+ - -/ M*a-p s (t)a. 



(20) 



Here 70 is the spontaneous emission rate given by 
70 = 4w 2 |d| 2 /3hc 3 , and er + , er_ are the standard rais- 
ing and lowering operators, respectively given by <r + = 
|1)(0|; cr_ = |0)(1| and 2N + 1 = cosh(2r)(2N th + 

1), M = -isinh(2r)e < *(2JV th + l), N th = (e^ -l)" 1 . 
Here Nth is Planck's distribution giving the number of 
thermal photons at the frequency w, r and <f> are squeez- 
ing parameters, d is the transition matrix elements of 
the dipole operator and c is the speed of light in vac- 
uum. We have neglected the Hamiltonian evolution part 
in eq™ ( 20 ) as we did in the evolution equation (JTJ . More- 
over, eq™(20) is in a Lindblad form and hence it corre- 



sponds to a completely positive map V [TTJ [12] . 

From equation) 20 1 we can find our V sq matrix which 
can be written as: 



a 








/9l 







ze' wt 








z*e iut 


ye iut 





/t 








v _ 



(21) 



where 

N{l + x 2 )+x 2 

2N + 1 
N(l~x 2 ) 

2N+1 
(N+l)(l~x 2 ) 

2N + 1 
jV(1 + x 2 ) + 1 

2iV+ 1 
exp[-7o(2iV+l)t] 

y = cosh f^y* J x 

( 7o a ^\ 
z = sinh ( I x exp |i$J 

a = sinh(2r) (2 N th + 1) 

The state corresponding to this channel V sq can be given 
by Choi-Jamiolkowski isomorphism [141 115) . The corre- 



ct- 

P 

v 

~2 
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sponding state is: 



{I®V sg )(\ct> + }(<t> + \)=M Sl 



a 













(3 


ze -iujt 










/' 





ye iut 








V 



(22) 



This matrix M S(J is a positive semi-definite matrix. The 
positivity under partial transposition, i.e, 



av- |z| 2 > 0, /3fi-y 2 > 0, 



equivalently 



27V + 1 



A(A+1) a f7o(2A+l)i 
27V + 1 Mn 1 



sinh' 



cosh" 



(23) 



>0, 
(24) 

> 0. 

(25) 



ensures separability. It is enough to show that equa- 
tion (25 1 is true in order to show that 



is separable as cosh' 



and sinh 2 (^f) < cosh" 2 



^ 7o(2iV+l)^ 



> 



the state 

2 



M. 



sinh z ( 
At t = this condi 



7 o(2JV+l)t 
2 



' 7o(2Af+lif 

,2 j 7q at 
^ 2 

tion ( |25[ ) is violated. At very large t, this condition is 
true if 2 A + 1 > a, i.e, 2 A + I > sinh(2r)(27V t/l + I). 
Now sinh(2r) < cosh(2r) which implies 2A + I > 
sinh(2r)(2A th + 1), since 2N + 1 = cosh(2r)(2A th + 1). 
This shows that there exists a finite time t at which M sq 
will turn from an entangled state to a separable state. 
This shows that one sided operation of V sq on \(j> + )((f) + \ 
evolves it into a separable state. Now using the factoriza- 
tion law for entangled states for d (g> 2 case (as was done 
we can say that all n-qubit states shows 



IV) 



m section 
ESD. 

When r = 0, i.e, when there is no squeezing then 
sinh(2r) = and 2A + 1 is always greater than 
sinh(2r)(2A t /j + l). This shows that when there is squeez- 
ing, the system takes more time to lose the entanglement. 
So instead of choosing the thermal bath one can choose 
squeezed thermal bath and delay the ESD process. 

At T = 0, ESD does not occur for a pure two-qubit en- 
tangled state if one considers the effect of V on one of the 
qubits where the initial state of the bath is a vacuum and 
the evolution is governed by equation ([I]). But for the case 
of squeezed bath, at T = 0, the mean occupation num- 
ber N is not zero and hence all the pure state show ESD. 
This implies that at zero temperature one can reduce the 
time of dissipation of the entanglement by switching on 
squeezing. Conversely, switching on squeezing at T =/= 
delays ESD. 



VI. QUBIT IN QUANTUM NON DEMOLITION 
(QND) INTERACTION WITH BATH 

In this section we will study the evolution of a sys- 
tem of qubits in QND interaction with bath. QND open 
quantum systems are those systems in which the sys- 
tem Hamiltonian commutes with the interaction Hamil- 
tonian. In this section we will be considering the evo- 
lution which Banerjee et al. has considered in [ID]. We 
can write the CP map Vqnd and the Mqjvd following 
equation (10) of reference [ID] : 



Vqnd — 



/ 1 \ 

e - iw * e ^(' iw ) 2 7(0 o 

e iu>t e -{h^) 2 ~f(t) q 

V o o o i / 



(26) 



M, 



QND 



1 e -*w* e -(fi w ) 7(*) \ 





V e ^* e -(M 2 7(*) o 1 



(27) 

where uj is the natural frequency of the system and j(t) 
is the time dependent spontaneous decay parameter (see 
[ID]). We can see that the matrix Mq N d is not posi- 
tive under partial transposition and hence under the one- 
sided action of Vqnd, |</ ,+ )(0 + | will not become separa- 
ble. Therefore this map will not show ESD. 



VII. SUFFICIENT CONDITION FOR ESD IN 
ANY FINITE DIMENSIONAL SYSTEM 

In this section we will derive the sufficient condition for 
ESD for any multipartite system where the dimension of 
each subsystem is d. For that purpose we will consider 
the state-channel duality for d ® d system. Consider a 
map $ acting on the B subsystem, where the state of the 
system AB is \x)ab — Si=i y/Pi\^)- The action of the 
channel on the state |x)ab can be written as: 

OU ® $) (Ix)^(xl) = 

A" : (® + \de (\x)ad(x\ ® PEC&)) \^ + )de ■ 
^E((I a ®$)(\x)ab(x\)) = 



E{M : {Ma ® 1 c )pac{M a ® Ic) 1 0- 



(28) 



where M A = Eft and p AC = (I ® $)(\<f> + )(4> + \). If 
$ turns out to be an entanglement breaking channel |16] 
(equivalently, if (IU ® $)(\4> + )ab(4 i+ \) is separable), then 
A" : (Ma ® ^c)pac(Ma ® ^cY ■ is also separable and 
hence, E((Ia ® $) (Ix)ab(xD) = 0- Thus separability of 
(I A ® $) (|0 + )as(V + I), for a given $, is sufficient for ESD. 

Now, based on the discussion in section |III[ we can 
claim that if we have a channel $ which shows ESD in 
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dtgid systems, then the same operation, acting locally on 
each subsystem, will show ESD in multipartite system as 
well. 



VIII. CONCLUSION 

In this paper we have shown that if the initial bath 
state is the thermal state, at finite temperature (T 7^ 0) 
all two-qubit states as well as multiqubit states show 
ESD. The same result has been shown for squeezed ther- 
mal initial states of the bath. We have found that under 
the one-sided action of the Markovian heat bath evolu- 
tion, one does not see any ESD at T — in two-qubit 
pure state where the initial state of the bath is the vac- 
uum state. This result is in sharp contrast with the case 
when the initial state of the bath is squeezed vacuum in 
which all two-qubit states show ESD. It has been shown 
that squeezing delays ESD in the finite temperature case 
but speeds it up at zero temperature. We have found that 
for the action of one-sided QND evolution with squeezed 
thermal state as the initial state of the bath, no two-qubit 



pure entangled state show ESD. Finally we have derived 
the sufficient conditions for the ESD for any multipartite 
system of identical (but distinguishable) particles. There 
are articles where one can find ESD in the systems evolv- 
ing under non-Markovian evolution |18) . 

After having the full understanding of ESD in Lind- 
blad type of evolution, the following problems become 
prominent: estimating the exact time to ESD for mixed 
entangled states; the factorization law for concurrence or 
any other useful measure of entanglement for two-qubit 
mixed entangled state; a general scheme for controlling 
ESD; a complete analysis of ESD in multipartite states 
and the quasi-covariance of dynamical evolution. One 
can have a full understanding of ESD only after all the 
problems have been solved. 
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